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Long-wave instability in shear flow 
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We prove here that the shortwaves in shear flows are always stable, whereas the shear flow 
instability is due to long-wave instability. For the unstable flows, an estimation of growth rate is 
obtained as toi < {p — l)y^/ir|(7max — f/min|, which is much better than the previous one obtained 
by H0iland and Howard. Comparing this with shortwave instabihty in Kelvin-Helmholtz (K-H) 
instability, it implies that the shear flow itself can damp the shortwaves without viscous effect. 
Further investigation shows that K-H instability is in essence of long-wave instability in shear flow. 
These results can be applied to verify numerical calculations. 

PACS numbers: 47.15.Fc, 47.20.Cq, 47.20.-k, 47.32. Cc 



Shear instability, caused by the velocity shear, is one of 
most important instabilities in the flows. Although the 
mechanism of shear instability is not full revealed yet, it 
is applied to explain the instability of mixing layer, jets 
in pipe, wakes behind cylinder, etc. Some simple models 
were employed to study the shear instability, including 
Kelvin-Helmholtz (K-H) model, piecewise linear velocity 
profile [H and continued velocity profiles 0, etc. To 
reveal the mechanism of shear instability, the stability of 
perturbation waves should be understood. In K-H model, 
the growth rate of the disturbance wave is proportion 
to the product of wavenumber and velocity shear, thus 
the shortwaves are more unstable than longwaves in K-H 
model. However, Rayleigh ^ found that the piecewise 
linear profile is linearly unstable only in a finite range of 
wavenumbers < k < kr, which means the shortwaves 
are stable in this case P, [^-Q. A contradiction d, Q 
emerges as shortwave and longwave dominate the shear 
instability in K-H model and piecewise linear profiles, 
respectively. This was explained either as viscous effect 
must be considered for shortwave Q or as longwaves do 
not "feel" the finite thickness of the layer [1, [J. 

To dispel the contradiction, we should find out which 
dominates the shear instability, shortwave or longwave? 
It is from the numerical simulation (see 0, S 13] and 
references therein) and some theoretical analysis (e.g. 
ToUmien Q and Lin Q) that longwave might be more 
unstable than shortwave. They also proved the long-wave 
instability subject to the velocity profile U is either sym- 
metric or monotone. We will investigate this problem 
following the way by Sun [Tol |. 

Moreover, a prior estimation of growth rate, which may 
be obtained from the investigation, is useful for unstable 
flows. For example, Howard's semicircle theorem Jll| has 
been used to validate the numerical calculations [3] . An- 
other useful result was found by H0iland [TJ and Howard 
[Tll | that the growth rate iVi must less equate to half of 
the maxim of vorticity, i.e., uJi < |?7'|max/2. But this 



estimation is too cursory to be used in applications. For 
example, U' is always great than zero even the veloc- 
ity profile has no inflection point. So this estimation is 
trivial for these cases. Here we will show a refinement es- 
timation of growth rate, which can be applied to general 
velocity profiles in parallel flows. 

The motivation of this short letter is to investigate 
these problems within the context of inviscid parallel 
flow. The aim here is to flnd out some general characters 
for unstable waves. 

To this purpose, long-wave instability in shear flows is 
investigated via Rayleigh's equation @-[al- For a paral- 
lel flow with mean velocity U (y) , where y is the cross- 
stream coordinate. The streamfunction of the distur- 
bance expands as series of waves (normal modes) with 
real wavenumber k and complex frequency uj = uJr + iuji, 
where uJi relates to the grow rate of the waves. The flow 
is unstable if and only if > 0. We study the stability of 
the disturbances by investigating the growth rate of the 
waves, this method is known as normal mode method. 
The amplitude of waves, namely cj), holds 
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where c = uj/k = Cr + ici is the complex phase speed. 
The real part of complex phase speed c,- = uJr/k is the 
wave phase speed. This equation is to be solved subject 
to homogeneous boundary conditions 



at y — a.b. 



(2) 



From Rayleigh's equation, we get the following equa- 
tions: 
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Rayleigh used only Eq.(|3]) to prove his theorem, i.e., a 
necessary condition for instability is U"{ys) = 0, where 
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Us is the inflection point and Us = U{ys) is the velocity 
at j/g. Before the further discussion, we need estimate the 

rate of WW^dy to /^^ M^dy, as Sun did 0. This is 
known as Poincare's problem: 



b pb 
'l|2. 



U'rdy = ti urdy 



(5) 



where the eigenvalue /i is positive definition for (j) ^ Q. 
The smallest eigenvalue value, namely //i, can be es- 
timated as /ii > (^3^)^- And an auxiliary function 

f{y) — u^u introduced, where /(y) is finite at 

inflection point. 

With the preparation above, we have such conse- 
quence. If —fiy) < Q < c», where Q is a positive con- 
stant, then the disturbances with shortwaves k > kc are 
always stable, where kc is a critical wavenumber subject 
to fc^ = Q — /ii. We will prove the consequence by two 
steps. At first, we prove proposition 1: if Cr = Us, the 
disturbances with shortwaves k > kc are always stable. 

Proof: Since U" = f{y){U-Us), Cr = Us and /(y) < 0, 
this yields to 
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> 0. 



This contradicts Eq.Q. So proposition 1 is proved. 

Then, we prove proposition 2: if Cr ^ Us , there must 
be = with k > kc- 

Proof: Otherwise if cf ^ 0, so according to Eq.Q and 
Eq.Q, for any arbitrary real number Ut which does not 
depend on y, it holds 

(ll0'lP + ^^ll'^lP) + ^^j|^^ll^lP]rf2/ = O. (8) 

But the above Eq.® can not be hold for some special 
Ut- For example, let Ut = 2c,. — Us, then there is {U — 
Us){U -Ut) < \\U-c\\^, and 
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For k > kc, this yields to 



\\m' + [k' + ^''^^_!^:h \\cbf}dy>o, (10) 



which also contradicts Eq.®. So proposition 2 is also 
proved. These two propositions are natural generaliza- 
tion of stabile criterion proved by Sun [Tsj . 

From the above two propositions, we can draw a con- 
clusion that the disturbances with shortwaves k > kc are 



always stable. This means that the shear instability in 
flows must be long-wave instability. Furthermore, the 
short-wave stability means that without any viscous ef- 
fect, the shortwaves can also be damped by shear flow. 
This mechanism is unlike the viscous mechanism that the 
viscosity has a damping effect on especially the short- 
waves. It implies that the shear flow itself can damp the 
shortwaves. This is very general and important conclu- 
sion, which explains why the instabilities found in shear 
flows are mostly long-wave instabilities. 

This result is also very important for numerical calcu- 
lation, which means shortwaves can be truncated in the 
calculations without changing the stability of shear flow. 
So the growth rates of longwaves have enough informa- 
tion for judging the stability of shear flow. On the other 
hand, the truncation of longwaves would probably change 
the instability of the shear flow. So the streamwise length 
scale must be longer enough to have longwaves for the nu- 
merical simulations in shear flows, such as plane parallel 
flow and pipe flow. Otherwise the instability of shear 
flow would be damped without long-wave perturbations. 

Then the growth rate of unstable waves can be esti- 
mated here by following the former investigation. Simi- 
lar to the assumption above, —fiy) < Q = p^/ii < oo, 
where 1 < p < 2, we will show that the growth rate is 

subject to UJi< {p- l)v7ir|C^max - C^minl- 

Proof: It is from Eq.® that gives 
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Substituting Eq.JS]) to Eq. pTj) and recalhng that < /i, 
this yields 



fk'Urdy< |^[^^|^^L_^-,^,]||0|Pdy. (12) 
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Multiplying above inequality ([T^ by cf , we get 
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u^t / Urdy< / h{y)Urdy 
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Suppose the maxim of h{y) is P^, then the growth rate 
is subject to 



This follows 



Hy)< 



, U"{U-Cr) , , 2 
\U-Cr)^+^, 



(15) 
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Substitution of /(y) into Ea. ([T5)) . gives 



,ip' -l)iU-Cr)'-C^ 
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When 

c2 = (p-l)([/-c.)^ (18) 

the right hand of Eq. ^T7\ get its largest value 

F2 = (p-l)Vi(t^-c.)^ (19) 
Then the growth rate must be subject to 

< (P- l)^/Ail|^/max - C^minl, (20) 

where J7min and J/max are minimum and maximum of 
U{y), respectively. And the wavenumber fcmax corre- 
sponding to the largest growth rate is 

fcmax = V(p - 1)^1. (21) 

So the result is proved. One should note that the fast 
growth rate uJi is only an approximation, but not a pre- 
cision one, so as to wavenumber fcmax. 

Comparing with the previous one by H0iland and 
Howard, the above estimation includes not only a more 
precision estimation about the the growth rate of unsta- 
ble waves but also the regime of unstable wavenumbers. 
This would be much helpful for validation in numerical 
calculations. 

As well known, the instability due to velocity shear is 
always associated to Kelvin-Helmholtz instability (K-H 
instability) , in which the disturbances of all wavelengths 
are amplified. According to K-H instability, the shorter 
the wavelength is, the faster the perturbation wave am- 
plifies. As the above investigation points out that the 
shortwaves are always more unstable than longwaves in 
this case, the contradiction at the beginning can be dis- 
pelled. 

An physical explanation 0, Q to the contradiction is 
that the K-H instability model has no intrinsic length 
scale, while Rayleigh's model has width of shear layer 



as length scale of waves. This can be noted from that 
Rayleigh's case reduces to the Kelvein-Helmholtz vortex 
sheet model in the long-wave limit fc < 1 HH, which 
is dramatically explained as longwave do not "feel" the 
finite thickness of the shear layer Q . Here we will show 
that this explanation can be extend to shear flows. It is 
from Eq. (|2ip that the fastest growing wavenumber /cmax 
is proportion to ^/JIi, then the corresponding wave length 
Amax is approximately 2(6 — a) / \/p — 1. So the thinner 
the shear layer is, the lager the fastest growing wavenum- 
ber is. It can be seen that the asymptotic case of infinite 
small shear layer leads to K-H instability. This is another 
evidence that K-H instability is in essence of long-wave 
instability. In this case, the K-H instability is an ap- 
proximation of shear instability when the wave length 
of perturbation is much longer than the width of shear 
layer. 

In summary, two general properties of shear instabil- 
ity are obtained in the investigations. Firstly, the short- 
waves are always stable in the continued profile flows and 
the shear instability is due to long-wave instability. This 
implies that the shear flow itself can damp the short- 
waves and that K-H instability is in essence of long- wave 
instability in shear flow. The result can be used to un- 
derstand the phenomena in hydrodynamics instabilities. 
Secondly, the growth rate of unstable flow can be es- 
timated as UJi < {p — l)y/JIi{Umax — Umin)- This esti- 
mation extend the previous result obtained by H0iland 
and Howard. Both results are important in numerical 
applications. The first one provides the estimation of 
unstable wavenumbers, and the second one provides the 
estimation of growth rate of unstable waves. These re- 
sults may be useful on both numerical calculation and 
stability analysis. 
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